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Abstract 

The decay amplitudes for — > DfD^ and — > D~^D~ have no factorizable 
contributions. We suggest that dominant contributions to the decay amphtudes arise 
from two chiral loop contributions and one soft gluon emission contribution. Then we 
determine branching ratios BR{W D+Dj) ~ 7 x 10"^ and BR(!B° D+D-) ~ 
1 X 10-3. 

Numerous experimental data coming from Ba Bar, Belle and Tevatron on B meson 
decays stimulate many studies of their decay mechanism. Through decades the factoriza- 
tion assumption has been used in calculations of the decay amplitudes. Recently, it has 
been shown [IJ that some classes of i?-meson decay amplitudes exhibit QCD factorization. 
This means that, up to a^/Tr (calculable), and kqcol'^h (not calculable), their amplitudes 
factorize into the product of two matrix elements of weak currents. Typically, the decay 
amplitudes which factorize in this sense are B ^ tttc and B — * Kit where the energy release 
is big compared to the light meson masses. However, for decays where the energy release is 
of order 1 GeV, QCD factorization is not expected to hold. Here we discuss the dominant 
contributions in B^ D^D~ and B^ D~^D~ j2]. At quark level these decays occur 
through the annihilation mechanism bs cc and hd —>■ cc, respectively (Fig.l). However, 
within the factorized limit the annihilation mechanism will give a zero amplitude due to 
current conservation, as for the —>■ K^K^ decay |B] . The axial part of the weak current 
might lead to non-zero factorized contributions if one of D-mesons in the final state is a 
vector meson D* . Such contributions are proportional to the numerically small Wilson 
coefficient Ci, which we will neglect in our analysis. In contrast, the typical factorized 
decay modes which proceed through the spectator mechanism, say B^ D^Dj , are pro- 
portional to the numerically larger Wilson coefficient C2. If one or both of charm mesons 
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in this decay are vector mesons, such amphtudes will give nonfactorizable chiral loop con- 
tributions to the process D^D^ due to i^'^-exchange. We determine these chiral 
loop contributions. 

There are also nonfactorizable contributions due to soft gluon emission. Such contri- 
butions can be calculated in terms of the (lowest dimension) gluon condensate within a 
recently developed Heavy Light Chiral Quark Model (HLxQM) J3], which is based on Heavy 
Quark Effective Theory (HQEFT) [5 . This model has been applied to processes with B- 
mesons in [HllZj. The gluon condensate contributions is also proportional to the favorable 
Wilson coefficient C2. We follow the standard approach |H] for non-leptonic decays where 
one constructs an effective Lagrangian Cw in terms of quark operators multiplied with 
Wilson coefficients containing all information of the short distance (SD) loop effects above 
a renormalization scale /i of order m^. Within Heavy Quark Effective Theory (HQEFT) 
the effective Lagrangian Cw can be evolved down to the scale /x ~ ~1 GeV P ITU]. 

The use of factorization is illustrated in the B^ —>■ D~^D~ decay: 

{D;D+\Cw\B')f = -(C2 + j^C,){D;\s^,^,c\0){D^\c^,b\W) , (1) 

The coefficients Ci^2 are Wilson coefficients for the operators containing the product of two 
left-handed currents. In our notation Ci = — ^KfeV^gflj, where the are dimensionless, 
Gp is the Fermi coupling, Vij are CKM parameters. Numerically, ai ~ 10~^ and 02 ~ 1 
at the scale fi = m^, and |ai| ~ 0.4 and \a2\ — 1.4 at /i ~ ~1 GeV CD]- Penguin 
operators may also contribute, but have rather small Wilson coefficients. 

The factorized amplitude for B^ — > D^D~ is presented in Fig. 1, and is given by 

{D-Dt\Cw\B')F = 4(Ci + ^C2){D;Dt\cZ^,CL\0){Q\TLYbL\B') . (2) 

Unless one or both of the Z)-mesons in the final state are vector mesons, this matrix element 
is zero due to current conservation: 

{DtD;\c^,c\0){0\dri,b\W) ~ fBiPD + PDT {DtD:\c^,c\0) = . (3) 

Our approach is based on the use of bosinized currents |2I and by using them we first 
write dawn the amplitude for 5° — > D^D^ . To calculate the chiral loop amplitudes we need 
the factorized amplitudes for B*^ DfD*^ and B^ D*^D*^, which proceed through 
the spectator mechanism as in Fig. |21 In this case the leading chiral coupling results from 
the coupling between a pseudoscalar meson H, vector meson H* a light pseudoscalar M 
(= TTjK,-)]), denoted by qa- After use of bosonized currents j2j, we obtain the following 
chiral loop amplitude for the process B^ —>■ D'^D^ from the Fig. 3: 

A{W DtD;y, = {v:jv:,) a(bI DjD;)^ ■ , (4) 

where the A{B'^ D^Dj)f stands for the factorized amplitude for the process B^ — >• 
D^Dj and the quantity is a sum of contributions from the left and right part of Fig. 
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Figure 1: Factorized contribution for D^D^ through the annihilation mechanism, 

which give zero contributions if both Df and are pseudoscalars. The double dashed 
lines represent heavy mesons, the double lines represent heavy quarks, and the single lines 
light quarks. 
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Figure 2: Factorized contribution for D^D^ through the spectator mechanism, 

which does not exist for decay mode B^ —>■ D^D~ we consider in this paper. 
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with uj = Mb/{2Md), a = [M|/(2M|,) - 1] and / being the vr decay constant. The 
function rix) is : 



r{x) 



X' 



== In + — , r(— x) = —r{x) + 



ITT 



(6) 



which means that the amplitude gets an imaginary part. Numerically, we find [2|: 

~ 0.12 - 0.26i . (7) 

The genuine nonfactorizable part for B^ D^D~ can, by means of Fierz transformations 
and identities for the product of two color matrices, be written in terms of colored currents 
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Figure 3: Nonfactorizable chiral loops for — > D^D^ . 
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Figure 4: Nonfactorisable contribution for — * DfD^ through the annihilation mecha- 
nism with additional soft gluon emission. The wavy lines represent soft gluons ending in 
vacuum to make gluon condensates. 



Within our approach, this amplitude is written in a quasi-factorized way in terms of matrix 
elements of colored currents: 

{DtD;\Cw\B')%F = 8C2 {DtD;\cL^,t''cL\G){G\dLi^t%L\B') , (9) 

where a G in the bra-kets symbolizes emission of one gluon (from each current) as visualized 
in Fig. 4. In order to calculate the matrix elements in (jHl), we have used j2] the Heavy Light 
Chiral Quark Model (HLxQM) recently developed in which incorporates emission of 
soft gluons modeled by a gluon condensate. Then we defined a quantity i?*^ for the gluon 
condensate amplitude analogously to in (jH) and (0) for chiral loops. Numerically, we 
determine j2] that the ratio between the two amplitudes is 

Rg ^ 0.055 + 0.162 , (10) 

which is of order one third of the chiral loop contribution in eq. ©. 

Adding the amplitudes R^ and Rq and multiplying with the Wilson coefficient [01 CSj 
a2 — 1.33 + 0.2i, we obtain the quantity: 

Rt = a2 {R^ + Rg) ^ 0.26 - 0.11? . (11) 

We have found that the amplitude for B^ is of order 15 — 20% of the factorizable 

amplitude for B^ D^D~, before the different CKM-factors are taken into account. 



Finally, we predict [2 that the branching ratios are 

BR{Wa D+D^) ~ 7 X 10"^ ; BR(fi° ^ D+D-) ^ 1 x lO'^ . (12) 

The current searches at BaBar and Belle might soon result in the limit on the rate 
— s> D^Dj. However, the 5° — > D^D^ mode will be accessable at Tevatron and later 
at LHC. 
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